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The conductance of disordered wires with symplectic symmetry is studied by the supersym- 
metric field theory. Special attention is focused on the case where the number of conducting 
channels is odd. Such a situation can be realized in metallic carbon nanotubes. The average 
dimensionless conductance (g) is obtained using Zirnbauer's super-Fourier analysis. It is shown 
that with increasing wire length, (g) 1 in the odd-channel case, while (g) — > in the ordinary 
even-channel case. It should be emphasized that the so-called Zirnbauer's zero mode, which 
has been believed to be unphysical, is essential for describing the anomalous behavior in the 
odd-channel case. 
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Phase-coherent electron transport in quasi-one- 
dimensional disordered wires is characterized by three 
universality classes. 1 The universality classes describe 
transport properties which are independent of the mi- 
croscopic details of disordered wires. The orthogonal 
class consists of systems having time-reversal symme- 
try without spin-orbit scattering, while the unitary class 
is characterized by the absence of time-reversal symme- 
try. The systems having time-reversal symmetry with 
strong spin-orbit scattering belong to the symplectic 
class. The conductance as well as its moment averaged 
over an ensemble crucially depends on which universality 
class it belongs to. For example, weak-localization (weak- 
antilocalization) effect arises in the orthogonal (symplec- 
tic) class, while such effects disappear in the unitary 
class. We can systematically analyze such differences by 
using two different nonpcrturbative approaches. 

The first is the scaling theory of Dorokhov, 2 and Mcllo, 
Pereyra and Kumar. 3 Let us consider an ensemble of dis- 
ordered wires of length L with N channels per spin. The 
conductance of a sample is expressed as G = (2e 2 /h)g 
with g = ^2 a= iT a , where {T a } denote the transmis- 
sion eigenvalues and the factor two in G corresponds 
to the electron spin. It is convenient to introduce A a = 
(1— T a )/T a . The statistical properties of g are determined 
by the probability distribution P(X\, A2, • • • , Ajv, L). The 
scaling theory is based on an evolution equation for P 
with increasing L, which is called the Dorokov-Mello- 
Pereyra-Kumar (DMPK) equation. The second is the 
supersymmetric field (SSF) theory of Efetov. 4 In this 
approach, we express the conductance (or its moment) 
in terms of electron Green functions and rewrite its en- 
semble average in an integral form over matrices Q con- 
taining commuting and anticommuting variables. Con- 
sequently, the problem is reduced to a field theory with 
the one-dimensional nonlinear a model. This approach 
enables us to calculate the average conductance and its 
variance by means of the super-Fourier analysis by Zirn- 
bauer. 5 ' 6 

The equivalence of the two approaches had been ques- 
tioned because, for the symplectic class, the Zirnbauer's 



result indicates (g) — ► 1/2 as L — ► oo while the DMPK 
equation gives (g) — ► 0. 7 It had been recognized that 
the nonzero limiting value of 1/2 is caused by the so- 
called Zirnbauer's zero mode, which appears in the super- 
Fourier analysis. This apparent controversy is resolved by 
Brouwer and Frahm. 8 They pointed out that the zero 
mode must be excluded since its parity is incompati- 
ble with an appropriate initial condition. They showed 
that the super-Fourier analysis gives us the result which 
agrees with the DMPK approach if we exclude unphys- 
ical modes, including the zero mode, and correctly take 
Kramers degeneracy into account. 

Recently, Ando and Suzuura 9 showed that metallic 
carbon nanotubes (CNs) have an anomalous transport 
property, which cannot be explained within the ordinary 
universality classes. They studied CNs of length L when 
the potential range of scatterers is larger than the lat- 
tice constant, in the case where the Fermi level lies in 
several conducting channels. Note that the number of 
conducting channels is odd in CNs. They proved that 
one perfectly conducting channel is present even when 
L — » 00. The present author 10 derived the DMPK equa- 
tion for this system based on the peculiar symmetry of 
the scattering matrix, 9 and found that the system be- 
longs to the symplectic class with an odd number of 
channels. A unit cell of a CN contains two sublattices, 
and this internal degree of freedom plays a role of a pseu- 
dospin. When the potential range of scatterers is larger 
than the lattice constant, there appears an additional 
symmetry, which drives the system into the symplectic 
class. 11 Thus, the symplectic symmetry in CNs has no 
connection with true spin as well as spin-orbit scatter- 
ing. Previously, the transport property in the symplectic 
class has been considered only in the even-channel case. 1 
However, the odd-channel case 9, 10, 12, 13 is very different 
from the ordinary even-channel case. 

In this letter, we consider the average conductance of 
disordered wires with symplectic symmetry by means of 
the SSF theory. Our attention is focused on the case 
where the number N of conducting channels is odd. We 
calculate the average conductance in the long-L limit by 
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using the super-Fourier analysis. In this approach, we 
expand the so-called heat kernel in terms of the eigen- 
functions of the Laplacian in the space in which the a 
model is defined. We point out that appropriate eigen- 
modes for the expansion differ between the even- and 
odd-channel cases. Consequently, unphysical modes in 
the even-channel case turn out to be appropriate in the 
odd-channel case. It is shown that the Zirnbauer's zero 
mode, which must be excluded in the even-channel case, 
becomes essential in describing the anomalous behavior 
in the odd-channel case. If we include the appropriate 
modes and correctly take Kramers degeneracy into ac- 
count, we obtain (g) ~ 1 + 2e _4L ^ in the odd-channel 
case when L » £, where £ = 2(N — 1)1 with the mean 
free path I. The result is in agreement with that obtained 
from the DMPK equation. 10 We hereafter treat only the 
symplectic class, and do not consider true spin. Our ar- 
gument closely follows that by Brouwer and Frahm. 8 
We first introduce a generating function F for conduc- 
tance, and show that the averaged F obtained by the 
SSF approach becomes equivalent to that by the DMPK 
approach. By inspecting this argument, we show which 
eigenmodes are appropriate in expanding the heat kernel 
and then obtain the average conductance. 

If the number of the conducting channels is even (i.e., 
N = 2m), the dimensionless conductance g is expressed 
as g = 2^^L 1 (1 + A a ) _1 , where the factor two represents 
Kramers degeneracy. Extending Rejaei's idea, 14 Brouwer 
and Frahm 8 introduced the generating function, 
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fo(M)=n/(M») 



(i) 



a=l 



with 



/(M) 



(1 + 2A + cos(0 3 + 04)) (1 + 2A + cos(fl 3 - fl 4 )) 
(1 + 2A + cosh 6»i f 

(2) 

Here, 9 — (81,63,84) and A = (Ai,A 2 ,--- , A m ). Once 
the ensemble average of F (8, A) is presented, we can 
calculate the average conductance, its variance and the 
density of transmission eigenvalues. 8 When the number 
of channels is odd (i.e., N — 2m+l), only one eigenvalue 
is not degenerate. We assign the Nth eigenvalue for it 
without loss of generality and obtain g=(l + Ajv) _1 + 
2^^ =1 (1 + A a )~ 1 . In this case, we consider the following 
generating function 



F(0,\) = f(0,\ N )* xFo(0,A), 



(3) 



instead of eq. (1). 

Before presenting the SSF treatment, we derive the 
evolution equation for the averaged generating function 
using the DMPK equation. 10 The ensemble average of F 
is given by 



(F(8, A))dmpk 



L 2 



dA r , 



dA 



N 



xS(\ N )F(8,\)P(\,L), (4) 
where P(X,L) = P(Ai,A 2 ,--- ,X m ,L). The probability 



distribution obeys the DMPK equation 
i(JV-l)Ap(A,L) 

m o q 



d\ a 



with the mean free path I and 

m m a—1 

J W = Il A cxIlIll A «- A ^ ( 6 ) 

c=l a=2 fe=l 

From eq. (5), we obtain the evolution equation, 

l(N - 1)^^(0, A))dmpk = (D X F(8, A)) D mpk, (7) 
where 

771 rj n 



d\ a 



In the limit of L —> 0, all the channels become perfectly 
conducting (i.e., A a = for any a). This ballistic initial 
condition results in 

r 1 vi a \\\ /cos6» 3 +cos0 4 \ 2m+1 (n , 

hm(F(8, A)) DM pk = [ 1 + CQsh9i ) ■ 0) 

In the even-channel case, the exponent 2m +1 is replaced 
by 2m. 

We next derive the evolution equation by using the 
SSF theory based on the model by Iida, Wcidcnmiiller 
and Zuk. 15,16 The derivation has been presented by 
Brouwer and Frahm, 8 so we only summarize it briefly. 
We assume that a quasi-one-dimensional wire consists 
of K small grains in series. Each grain has 2M electron 
states including Kramers degeneracy The wire is coupled 
to two ideal leads. We assume that the left (right) lead 
contains N\ (N2) conducting channels. The Hamiltonian 
H of our system is described by the matrix elements 
H l j s ^ VSi/ , where i, j denote the grains (1 < i, j < K), /x, v 
denote the electron states (1 < fi, v < M) and s M , s u 
denote the pseudospins. The elements H" VSi/ , which 
describe the electron states in the zth grain, are indepen- 
dent Gaussian random variables with zero mean values 
and 

2 

^ / 

- W^'( i(7 !/)v s f( i0 '!') s -^')' ^ 10 ^ 

where a y is the Pauli matrix. The coupling between the 
adjacent grains are described by H l J s ^ us ^ (j — i ± 1), 
which are another set of Gaussian random variables with 
zero mean values and 



TTl] TTjl 



(11) 



H^s li us t ,H^, s ^ I/ , s ^ / — ^-^5pp'5 ulJ i{i(jy) s ^, tSfji (i(jy) Sv , Sil . 

(12) 

All other matrix elements H Z J S vs ^ vanish. The coupling 
to the leads is described by a fixed 2KM x (N1 + N2) ma- 
trix W. Its elements W^ vo , where a denotes the modes in 
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the leads, are nonzero only for i = 1 and 1 < a < N± , and 
i~K and N\ + 1 < a < Ni+N 2 - We assume that the or- 
thogonality relation tt^^ W^W^j = wi^ q6 holds 
for i = l,K. 15 We decompose W &sW = W1+W2, where 
W\ and W2 contain the block with {W* a a } and that 
with {W^ „}, respectively. Brouwer and Frahm proved 
that the generating function eq. (3) can be expressed in 
the form, 

F(§, A) = Sdct^s ^ F i 8 _ Hl 8 + ittWxWIq 



has the same symmetry as the heat kernel, we eventually 
find that (F(§, A)) SSF obeys 8 

8(v 2 /vi) 2 ^))ssf = A S {F(6, A)) SSF 

with the initial condition 



lim(F(0,A))ssF= / dQ'f 1 (Q,Q')f2(A,Q'). 

s— »0 



(21) 
(22) 



mW 2 W$A), 



The Laplacian is defined as 



(13) 



where Sdet stands for the superdeterminant of a super- 
matrix and E-p the Fermi energy. The supermatrix Q is 
expressed as Q = T~ 1 AT with A = diag(l 4 , — 14) and 



.3 = 1,3,4 

with 

J(0) = sin6» 3 sin6»4sinh 3 6ii 
1 



(23) 







cxp 















u 
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- (0! +i Sl 6> 3 + is 2 4 ) 



(24) 



where u and w are 4x4 unitary supermatrices and 



(14) 
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( 01 
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(15) 



x ] J sinh 

Si ,S2 = ±1 

We prove that the SSF approach is equivalent to the 
DMPK approach. It is convenient to introduce D^°\ 
which appears in the evolution equation for the even- 
channel case, 8 given by 



d\ a 



(25) 



We can rewrite F in terms of a Gaussian integral over 
supervector ip as 



F(§,\) = J 



with J (A) = n^L2ll6=i l^o - A 6 | 4 . Using the rela- 
tion Df ] F (e, A) = A § F (9, A), derived by Brouwer and 
^A(jS F l 8 - Hl 8 + mW-iWlQ Frahm, we obtain 



+ i7rWW 2 T A + ieA U 



(16) 



where e is a positive infinitesimal. After performing the 
standard manipulations, 4,6 ' 15 we find that in the limit 
M — > 00, the averaged generating function is given by 



(F(6,X)) 



SSF 



where 

W(Q U Q K ) = 



dQi 



K-l 



dQ K fi(Q,Qi)f2(A,Q K ) 
x W(Qi,Q K ), (17) 



A § F(6,X) - [Df +J28g(0,\a)j F0,X) 
(A e f(6,\ N )t)F (6,\), 

dlnf(§,X N )i dlnf0,X) 



+ 



where 
6g0,X) 



.3=1,3,4 



09, 



(26) 



(27) 



We can show that A § f(9, X n ) 1 / 2 \ Xn=q = and 



Y[dQj 



i=2 



D^+^9(0,\a)j F(§,X)\ Xn=0 = D- x F(9,X)\ Xn=q , 

(28) 



if-i 



x exp -(w 2 /fi) 2 E Str(Q fe Q fc+ i) 



k=i 



h(Q- ^ ) = <'xp( -^Strln(l + xQQ,) 



h{Q,Qx) =exp 



AT 2 



Strln(l + xQQ K ) 



(18) 

(19) 
(20) 



\ 0=1 / 

which result in 

A § F&X)\ Xn=q =D x F(9/X) 
Using this relation, we obtain 

l(N- A)) DMPK 



=0' 



(29) 



Here, Str denotes the supertrace for a supermatrix. As- 
suming that v\ -C v 2 , we replace the sum in cq. (18) 
with an integral. 6 ' 15 Accordingly, the discrete number k 
is reduced to the continuous variable s. The propagator 
W(Q,Q'), called the heat kernel, obeys the heat equa- 
tion in the supersymmetric space. Noting that f\{Q, Q') 



„ = 



„ = 



Y[dX 

a=l 
m 

n ,lv 



dX N 6(X N )P(X, L)D k F(9, A) 



0=1 



dX N 6(X N )P(X, L)A S F(§, A) 



A§(F(9, A))dmpk- 



(30) 
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Thus, the evolution equation for (F(9, A))dmpk is equiv- 
alent to that for (F(0, A))ssf if we set 

/ [ \ ii ■ 

(31) 



l(N-l) 



2L' 



Now we focus on the initial condition. We assume that 
the coupling to the two leads is ideal, which is achieved 
by setting x = 1. Brouwer and Frahm showed that the 
ballistic initial condition for the DMPK approach (i.e., 
A a = for any a) can be reproduced if we fix Ni and 
take the limit AT 2 — > 00. We trace the same steps. The 
limit N 2 -*■ 00 results in / 2 (A, Q') -> S(Q',A). Thus, the 
right-hand side of cq. (22) is reduced to /i(Q,A). We 
finally obtain 

lim(F(0, A)) SSF - cxp f-^Strln(l + QA) 



COS 63 + COS 64 

1 + cosh Q\ 



N,. 



(32) 



This is equivalent to eq. (9) if we set Ni = N = 2m + 1. 
Thus, the SSF approach is equivalent to the DMPK ap- 
proach if we accept the treatments described above. This 
indicates that the average conductance as well as its vari- 
ance is the same whether computed from the SSF ap- 
proach or from the DMPK approach. It should be em- 
phasized that eq. (32) becomes the initial condition for 
the even-channel case if we set N\ = 2m. 

We here briefly summarize the Zirnbauer's super- 
Fourier analysis 5,6 for the thick- wire limit which is de- 
fined by N —> 00, L/l — > 00 with a fixed Nl/L. In this 
analysis, we expand the heat kernel W(Q, Q') in terms of 
the eigenfunctions of A^. The eigenvalue equation is 5,6 ' 8 



A, 



(0) 



1 



£(A,rai,n 2 )^A,ni,7i2(0)- ( 33 ) 



The eigenvalues are e(X, m, n 2 ) = A 2 



lo — 1 , where 



rii, n 2 = 1, 3, 5, • • • and A is a positive real number. In 
addition, there arises the zero mode e(i, 1,1) = 0, and 
the subsidiary series e(i, n, n — 2) = 2(n — l) 2 , where 
n = 3, 5, 7, • • • . Using the resulting expansion, Zirnbauer 
obtained the expression of the average conductance, 

(#) Zirnbauer = T(i, 1, l)/2 

+ {T{i,n,n-2)+T{i,n-2,n))/2 

n=3,5,7,- 

+ 2 4 V / dAA(A 2 + 1) tanh(7rA/2) 

m,n2=l,3,5," ^° 

x nin 2 (A 2 + n\ + n\ - 1) 



x Q (-1 + itrA + ami + ami) 1 

(T,(Tl,0 - 2=±l 



x T(A,m,n 2 ) 



(34) 



with T(A,m,n 2 ) = c -( A2 + n ?+»l-i)(i/2«). This indicates 
that (g) — *• 1/2 when L — > cxo. 5,6 This apparently con- 
tradicts the DMPK approach, which gives (5) — > in 
the even-channel case 1 and (g) — ► 1 in the odd-channel 



case. 
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Note that the difference between the even- and odd- 
channel cases is not reflected in the Zirnbauer's argu- 
ment. We show that the even-odd difference can be taken 
into account through the initial condition for the heat 
kernel. Let us consider an operation P defined by 

P<j>{6l, 3 , 4 ) = 0(01, 03 + 7T, 04 + 7T), (35) 

which is introduced by Brouwer and Frahm. 8 The Lapla- 
cian Ag commutes with P, and the eigenfunctions have 
either even or odd parity. The even (odd) parity means 
that the eigenvalue of P is +1 (—1). Zirnbauer considered 
the limit Ni — 7V 2 — > oo. If the same limit is taken in our 
argument, we find that (F)ssf = W(Q,A) holds 8 since 
fi(Q,Q') is reduced to the delta function S(Q,Q') when 
Ni — > oo. This indicates that we can obtain the heat ker- 
nel by expanding (-F)ssf in terms of the eigenfunctions of 
Ag and then taking the limit Ni — > oo. By employing this 
method, we can obtain W(Q, A) which reflects the even- 
odd difference. The crucial point is that before taking 
the limit Ni — > oo, the parity of eq. (32) is well defined; 
its parity is odd (even) when Ni = 2m + 1 (Ni = 2m). 
Accordingly, we must select only the eigenmodes hav- 
ing odd parity in the odd-channel case, while the even- 
parity modes are appropriate in the even-channel case. 
However, once the limit of Ni —> oo is taken, eq. (32) is 
reduced to 



\im(F )ssf 

s— >0 



S(Q,A), 



(36) 



which is equivalent to the Zirnbauer's initial condition, 
and we cannot derive such a parity-selection rule from 
it. Thus, the even-odd difference disappears. In order to 
correctly calculate the conductance, we must take care 
of this parity-selection rule. Brouwer and Frahm pointed 
out that the odd-parity modes are unphysical in the even- 
channel case. However, they were not aware that the odd- 
parity modes become essential in the odd-channel case. 

In the even-channel case, the above argument tells us 
that we must neglect the eigenfunctions with odd par- 
ity. The zero mode as well as its subsidiary modes has 
odd parity, 8 so they must be excluded. By omitting all 
the odd-parity modes and correctly taking account of 
Kramers degeneracy, Brouwer and Frahm obtained 



W/ 9 \2L 



e 2 «. 



(37) 



which agrees with the DMPK approach. 1 In contrast, we 
must exclude the even-parity modes in the odd-channel 
case. This is the crucial difference between the ordinary 
even-channel case and our odd-channel case. If we omit 
all the even-parity modes and take Kramers degeneracy 
into account, we obtain 



<0>=T(i,l,l) + 



(T(i,n,n-2)+T(i,n-2,n)) 



n=3,5,7,- 



+ 2 5 



ni,n 2 

711+712 



-1,3,5,- Jo 
E0(mod4) 



dAA(A 2 + 1) tanh(7rA/2) 



x n!n 2 (A 2 +n\ + nl-l) 

x || (-1 + icrA + ami + <7 2 n 2 )" 

<T,(T1 ,(72 — ±1 
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xT(A,m,n 2 ). (38) 

When L — ► oo, the leading contribution to (g) arises 
from the zero mode, which yields (g) — > 1. This indicates 
the presence of one perfectly conducting channel. The 
deviation (Sg) = (g) — 1 from the quantized value is 

(Sg) ~2e~ ± t, (39) 

which comes from the subsidiary modes with n = 3. This 
clearly indicates that the zero mode as well as the sub- 
sidiary modes is essential in describing the anomalous 
transport property in the odd-channel case. The char- 
acteristic length scale for (Sg) is £ dd = £/4, while that 
in the even-channel case is obtained as ( cvcn = 2£ from 
eq. (37). Thus, £ dd is a factor 8 shorter than £ eve n- This 
result is in quantitative agreement with the DMPK ap- 
proach. 10 
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